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In this paper, a description, and analysis of a novel 3-D dimension 
hyperchaotic system is implemented. The proposed system oscillation is two- 
order autonomous and consisted of a nine-term and symmetric oscillation w.r.t 
x-axis. It is proved analysis by Kaplan-York dimension, waveform analysis, 
phase portrait, and Lyapunov exponent. This work-study stability and 
equilibrium point and Routh stability criteria produced that the new system 
has one unstable point from the type saddle-focus point. One of the 
characteristics of the proposed system is hyperchaotic since this system has 
two Lyapunov large than zero. This system is applied to generate a chaotic 
matrix;¢.1¢6 (S-box) based in advanced encryption standard (AES) algorithm 
for text encryption and gives a high level of security. In addition to the 
description, and analysis S-box. Therefore. the proposed algorithm is satisfied 
the high randomness of entropy value and passes the National Institute of 
Standards and Technology (NIST) parameters and another test. Mathematica 
and MATLAB programs simulated some results. 
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1. INTRODUCTION 

In the past few decades, researchers have studied chaotic systems, chaos theory, and applications 
based on chaos a focal study point on the 3-D chaotic systems since the discovery of the first chaotic attractor 
by Lorenz in 1963 [1]. In recent years several chaotic systems have been established including Chen and Ueta 
[2], Zghair et al. [3], Hajipour and Aminabadi [4], ElBeltagy et al. [5], and lu systems [6], also self-excited 
and hidden attractors are two recent categories of chaotic attractors, attraction has a basin that does not cross 
any tiny neighborhoods of an equilibrium point which is called a hidden attractor, in contrast to [7], an attraction 
basin which has crossed tiny neighborhoods of an equilibrium point which is called a self-excited attractor. 
Zhang et al. [8] worked on a novel 3D chaotic system by applying a nonlinear cross-product to the 2"? equation 
of the Rucklidge system they found that the worked system has extremely complex dynamical behaviors [8]. 
Lassoued and Boubaker [9] provides a general survey about newly designed chaotic and hyperchaotic systems, 
this involves an analysis of the dynamical behaviors of these complex systems Salih et al. [10] designed a new 
approach by using a 3D Logistic map as a proposed key through the replacement of the predetermined XOR 
operation in advanced encryption standard (AES) algorithms. Faghani et al. [11] designed a new 3D chaotic 
system depending on the basic formula of quadratic Jerk systems having terms, and other systems having a 
stable node equilibrium point. Patra and Banerjee [12] prove that 3D systems have attractors such as higher 
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dimension torus and hyperchaotic which cannot be seen in 2D. Abdullah and Abrahim [13] try to improve the 
Rivest-Shamir—Adleman (RSA) algorithms by applying a 3D chaotic system with an experiment characteristic 
designed for this purpose. The new proposed system analysis and describes complex dynamical behaviors by 
proving the Lyapunov exponent spectrum, phase portraits, symmetric, dissipation, entropy, SDIC, and Routh 
stability criteria, and waveform analysis of the proposed system. This article introduced a novel suggestion for 
improving the security of AES as a chaotic secure communication algorithm. The improvement is driven by 
generating an XOR matrix 16x16 S-box table form for the AES algorithm using the proposed system. This 
article's body is organized as follows: section 1, analysis and describes the mathematical modeling of chaos 
behavior, the stability and equilibrium of a new 3D proposed system, and introduced Kaplan- Yorke dimension, 
Lyapunov exponent and waveform test. Section 2 describe proposed metode. Section 3, the proposed system 
generated a chaotic XOR matrix ,¢,.4¢. Section 4, security analysis Entropy and National Institute of Standards 
and Technology (NIST) test result and discussion. The conclusion is presented in the final section. 


2. THE PROPOSED METHOD 

In this part overview analysis and describes a new 3D proposed system behavior. Also, stability and 
equilibrium, Kaplan- Yorke dimension, Lyapunov exponent, and waveform test of a new 3D proposed system. 
Addition generated S-box based on 3D proposed system; i) analysis and describes new mathematical modeling 
of chaotic behavior and ii) the proposed system generated a chaotic XOR matricy¢ x16. 


2.1. Analysis and describes mathematical modeling of chaotic behavior 

Considering the 3D hyperchaotic system the proposed system is two-order autonomous and consists 
of nine terms where xy, xz are the quadratic cross-product nonlinear term and it has two-term, one of them. 
sine and another cosine function, the proposed hyperchaotic 3D system is written as expressed: 


= = —n x — 0z? + B Cos(y) 
Z = pSin(x) — y — xz (1) 
dz 
a Oy t UXY — Nz 

The proposed system has seven positive parameters f, 0,7, 2, p , w and p, and three initial conditions 
(x(0) y(0), z(0) ), where (x,y,z)" € R°, a strange attractor and hyperchaotic behaviors generated by the 
proposed system when (,0,n,2,p,@, W=(1.5, 1.38, 11, 2.5, 30, 15, 5) and (x(0) y(O) 2(0)),=(3.6,1,2.6), 
therefore the keyspace for this system is very high and by the coordinate transformation: 
(x, Y, Z) — (x, —y, —Z), the proposed 3D hyperchaotic system (1) is invariant. Therefore, 3D hyperchaotic 
system (1) is a symmetric oscillation w.r.t x-axis. Figure 1(a) phase partial for the proposed system in 2-D plan 
space (x,y), Figure 1(b) phase partial for the proposed system in 2-D plan space (x,z), Figure 1(c) phase partial 
for the proposed system in 2-D plan space (y,z), and Figure 1(d) attractors simulated in mathematica which 
display phase partial for the proposed system in 3-D (x, y, Z) space. 


2.1.1. Stability and equilibrium of a new 3D proposed system 
By setting x = y = z =0, to solve the system (1): 


0 = =y x — 0z? + B Cos(y) 
0 = p Sin(x) — y — xz 
0 = wy + xy — Nz 


We obtain by substitution (6,0,7, 2, p, w, w)=(1.5, 1.38, 11, 2.5, 30, 15, 5) one equilibrium point v 
(0.00012,0.0037,0.0223), the linearized for the proposed system (1) at v with Jacobian matrix is: 


=n —Sin(y) —2z0 —11  -0.00367 —0.0615 
j = |—z + pCos(x) —1 =x |=|29.9776 —1 —0.00012] , let |[AZ — j| = 0 
yu xuto =R 0.0555 5.0018 —2.5 


The eigenvalue of v is gained as 1,=-3.34 — 6.91 i, A, = —3.34 + 6.91i and A3 = 14.188. The 
system at the point v is a saddle point, and the saddle-focus v has two dimensions stable manifold and a one- 
dimension unstable manifold, so this equilibrium point is unstable. Therefore, the characteristic equation is 
—2? + (-1-n-2)/? (—n — 2yzuo + x(—xp — w) — 2 — nO — zBSin(y) + BpCos[x]Sin(y))A + 
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ZBOQSin(y) + BpNCos(x)Sin(y) = 0 which implies that: 


x(t) 


x(t) 


yE) 


(c) (d) 


Figure 1. Attractors portrait (a) (x, y), (b) (z, x), (c) (x, z) in 2D space, and (d) (x, y, x) in 3D space 


(2) 


372.971 — 55.9084 — 14.5 2? — A?=0 
=, 


We obtain the Routh—criterion [14] from (2), aọ = 372.971, a, = —55.908, a, = —14.5 and, a3 
— 8% — -30.1859 hence the proposed system (1) is unstable also we conclude the Hurwitz — 


oa ^ 

criterion by formed A1= @ = -14.5 < 0, A= is i) = 0.001118 > 0, 
az aA 0 

A= |a} a, 0] =0.000044 > 0, hence the proposed system (1) is unstable. 
0 ay a 


2.1.2. Kaplan-Yorke dimension, Lyapunov exponent and waveform test 
Its estimated divergence and convergence of neighboring paths are determined by nonlinear chaos in 


a dynamical system called the Lyapunov exponent also defined [3], [15]. Therefore by the three Lyapunov 
exponents when (p,o,n, Q, p,w, u ) = (1.5, 1.38, 11, 2.5,30,15,5 ) and (x(0) y(0) z(0)),=(3.6,1,2.6), are 
calculated by theoretical and numerical Alan Wolf analysis with step size 0.05 used Jacobian matrix after 
10000 iterations, hence L}, = 2.74 > 0,L,, = 0.0335 > 0, and L}, = —17.249 < 0, since an Ly, >0 
and yo, Ly, < 0 after the order Ly, > La, > Ly, , we obtain the proposed system (1) is hyperchaotic, and an 
upper bound for fractal dimensions and Hausdorff which is called is Kaplan-Yorke dimension [16] i.e. 
KY,, = j+ m= Ly, hence as a result of this, KY,, = 2.1608. Figure 2 shows the Lyapunov diagram 
Aj+1 

have Ly,, > 0,L,, <0. 
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Waveform test: the waveforms [15] for the 3D proposed system (1) are included in Figure 3 in the 
time domain, as a result of the oscillation's aperiodicity, their time-domain features are noncyclical. 
Figure 3(a) shows time vs x in the proposed system (1). Figure 3(b) shows time vs y in the proposed system 
(1). Figure 3(c) shows time vs x in the proposed system (1). 


Dynamics of Lyapunov exponents 


Lyapunov exponents 


oO 100 200 300 400 500 600 700 800 900 1000 
Time 


Figure 2. Lyapunov diagram by setting ( 8,0,7, 2, p, œw, u )= (1.5, 1.38, 11, 2.5, 30, 15, 5) 


time aginest x(t) time aginest y(t) time aginest z(t) 


(b) (c) 


Figure 3. Waveforms for new system w.r.t (a) x, (b) y, and (c) z against time 


2.2. The proposed system generated a chaotic XOR matrix 4¢,1¢ (s-box) 

Cryptography science comes to play an integral role in many applications during the ancient time in 
many wars. It is used mainly in the military to send a secret message from one side (sender) to another side 
(receiver) [17], [18]. The sender hides the original data for protecting it from unauthorized access. AES 
algorithm is widely used these days as a specification for the encryption by U.S NIST. It considers a much 
stronger than the other encryption algorithms such as DES and DES because of its block cipher, key size with 
128/192/and 256 bits. The block cipher for the encrypted data is 128 bits. The basic definition of encryption is 
a method for protecting information such as text, documents, images, and files from changing, modification, 
and eavesdropping other than the intended recipient. The secret key is used by authorized people to encrypt 
and decrypt the data, where the sender uses one type of mathematical formula to convert the plain text to cipher 
text in a format that cannot be understood by the jeopardized users. The receiver performs the reverse process, 
taking the ciphertext and converting it to the original using the key agreed upon between the two parties. 
Currently, cryptography plays a major role in securely transmitting and receiving messages in online banking, 
telecommunications, the IoT, fog computing, and cloud computing [19]. The encryption process is classified 
into two types: the first one is using public and private keys, and the second is using the shared key between 
the two parties. 
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2.2.1. Proposed algorithms: the proposed XOR matrix,,¢,1¢ consist of several points 

Input the parameters B, 0, 9, Q, p, w and y . and (x(0),y(0),z(0)) which represent the initial conditions 
to generate three sequences {( x;,y; and zi); i = 1, 2,..., n} according to the proposed system (1). XOR 
between ( x; , y; and zi) after their conversion to unsigned 16-bit integer numbers by following the equation: 


x_i = uint16(x_i * 10*10,256) 
y_i = uint16 (y_i * 10*10,256) 
z_i = uint16 (z_i 10*10,256) 


generated XOR. 

The S-box new has the following properties: i) bijective, ii) the number between (0 — 255 ) without 
periodic in every row and every column, and iii) 127.5 is equal to the average of all matriXx16»16 entries. 
Figure 4 show the overall proposed system for the dynamic behavior with AES encryption. The only difference 
between a proposed XOR matrix (16*16) and XOR matrix (16*16) in AES [10], [20] is that the proposed XOR 
matrix (16*16) uses integers produced from a new 3D hyperchaotic system (1). Table 1 shows the XOR 
matrix (16*16) S-box generated based on the new 3D hyperchaotic system. 


Data Input in 
any Format 


4 


z AES S-box generation based chaotic dynamical 


s Input the parameters 


B,0,n,2,p,w and p .and 
the initial conditions (x(0), 
: " y(0),z(0)) 


Initial Key Ki 


Apply 3D chaotic system (1) 
For a chaotic keys 


y= uint16(y;*10^10,256) 
z= uint16(z;*10^10,256) 
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Figure 4. The overall proposed system for the dynamic behavior with AES encryption 
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Table 1. The proposed XOR matrix,,¢,1¢ S-box generation 
0 1 2 3 4 5 6 vi 8 9 10 11 12 13 14 15 
163 248 107 155 19 223 188 133 36 150 100 245 97 137 113 193 
124 25 141 86 94 18 234 229 143 87 147 142 168 242 177 3 
70 115 152 157 179 2 215 60 30 191 238 197 101 35 151 186 
131 232 130 146 84 253 74 126 182 198 254 68 110 148 240 57 
190 203 196 136 73 62 169 237 56 14 158 88 139 217 
105 189 5 216 181 171 129 108 6 76 221 20 128 47 132? -233 
120 52 246 219 75 206 95 67 99 69 125 61 78 64 211 104 
83 183 145 13 184 26 22 252 194 72 121 162 180 15 50 144 
10 11 43 27 98 202 247 40 138 226 170 106 38 59 244 149 
173 24 167 122 160 161 134 90 116 103 153 187 66 41 37 71 
10 227 «218 235 23 119 117 214 208 32 77 212 135 0 45 192 55 
11 46 49 249 7 9 53 82 109 236 228 39 255 204 118 200 63 
12 16 102 209 172 44 85 241 79 12 96 165 34 159 175 166 220 
13 251 28 174 33 127 65 222 58 243 17 231 250 195 42 54 4 
14 140 80 164 225 93 213% +1 29 89 239 81l 185 156 111 112 31 
15 210 224 201 207 51 91 21 92 48 199 205 8 176 178 230 114 


CMONIDNWKRWNHEO 
= 
Ur 
p 
= 
N 
oO 


3. THE PERFORMANCE ANALYSIS OF S-BOX 

The S-box generator is evaluated using many criteria. The nonlinearity provides S-box with 
uncertainty output, which provides counter measure to differential and linearity cryptanalysis assaults, while 
the strict avalanche criterion (SAC) is used for explaining the probability of each bit changing when one bit in 
a Boolean functions input changes is 0.5. Time is calculated to show the delay for the generation of the proposed 
matrix_(16*16) (S-box) based on a new 3D hyperchaotic system, In this work a new 3D hyperchaotic proposed 
system using (8,0,n, 2, p, w, u )=(1.5, 1.38, 11, 2.5, 30, 15, 5)and (x(0) yO) 2(0)),=(3.6,1,2.6), therefore the 
keyspace [21] for this system by taking 15-digit precision is more than 101°° ~ 24°, it is very high. Also, the 
proposed S-box was measured and examined concerning the indicated features: 


3.1. Nonlinearity 

The nonlinearity is shown in [22]. We observe that the proposed matrix;¢.;¢ (S-box) nonlinearity 
rates are 110, 107, 108, 110, 108, 109, 110, and 109. We found that miny; = 107 and average y; = 108.785, 
all component functions had very large nonlinearity value. it is obvious that, the proposed (S-Box) matrix,¢.16 
has excellent nonlinearity behavior. 


3.2. SAC 

The SAC which is shown in [22] also tested in this work for a matrix16+165-box generation based on 
a new 3D hyperchaotic system. Table 2 shows the main value of this test. It shows that the rate value for 
SAC matrix1¢.1¢6 S — box = 0.5016, therefore any change in the input then the matrix ;¢,.1¢ (S-box) has good 
SAC and a decent avalanche. Table 3 shows the comparative study of the proposed matrix,.¢,;¢ (S-box) with 
other references w.r.t rate value of SAC. 


Table 2. SAC for the proposed XOR matrix16+16 S-box generation 
4 5 6 


NO. 1 2 3 7 8 
1 0.58 0.5 0.45 0.5 0.48 0.47 0.58 0.52 
2 0.44 0.48 0.48 0.47 0.5 0.5 0.48 0.53 
3 0.48 0.56 0.48 0.56 0.5 0.44 0.56 0.52 
4 0.5 0.52 0.48 0.5 0.48 0.52 0.47 0.52 
3 0.47 0.45 0.5 0.55 0.5 0.56 0.52 0.47 
6 0.44 0.48 0.47 0.53 0.47 0.47 0.52 0.52 
7 0.58 0.47 0.55 0.42 0.55 0.58 0.45 0.53 
8 0.48 0.53 0.53 0.45 0.48 0.44 0.56 0.53 
Table 3. Comparative study of the S-box w.r.t rate value of SAC 
S-box Proposed S-box Ref [22] Ref [23] 
rate value for SAC 0.5016 0.4976 0.4981 


3.3. Time performance 
One of the most components of any security is calculated to speed time. Therefore, we computed the 
time for the generation of the proposed matrix16»16 (S-box) based on a new 3D hyperchaotic system it was 
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discovered the technique on average a very small-time rate=0.275. The hardware used to implement this work 
is Intel Core 17, window 10, RAM 8 GB, 64-bit MATLAB 2019b, and mathematica 11. 


4. RESULTS AND DISCUSSION 
The main measurement to evaluate this work is entropy which is defined by Shanon law [24] in the 
(3): entropy is implemented in most cryptography systems to evaluate the randomness of the encrypted text. 


H(x) = —Y%2 P(xi)log,(p(xi)) 372.971 — 55.908 — 14.5 42 — 23 =0 (2) 


Table 4 shows the main entropy value comparison between AES and chaotic dynamical behavior AES 
(CDB-AES). It is clearly shown from the table the entropy of the proposed method is close to the entropy of 
the AES algorithm. However, the ideal entropy value ought to be eight. If the entropy of the cipher shows a 
result less than eight, there exists a percentage of obviousness and the security will be a violation. The result 
of entropy is calculated on different message bits of the proposed and the AES algorithms using 256 different 
chaotic keys. Table 4 shows the comparison between the value entropy for both the proposed and the AES 
algorithms. The average entropy of the proposed algorithm (CBD-AES) is 7.746414. While the standard AES 
is 7.731757. 


Table 4. Comparison between AES and improved CDB-AES 


#No AES/ entropy CDB-AES/entropy Message Length/bits Chaotic K-length/bits 
1. 7.2806 7.6525 128 256 
2. 7.9077 7.7990 256 256 
3. 7.9600 7.9866 512 256 
4. 7.9009 7.9326 1024 256 
5. 7.8090 7.5693 2048 256 
6. 7.6796 7.6994 4096 256 
7. 7.5845 7.5855 8192 256 


4.1. National Institute of Standards and Technology test 

The randomization of new 3D hyperchaotic keys for the proposed XOR matriXx16+16 the process is 
tested by the national institute of standards and technology package for binary bits sequences [25]. Table 5 
results show that the new system is passed successfully in all parameters. Therefore, the proposed XOR 
matrix,¢+16 has produced a sequence with excellent, arbitrary-satisfying features in all different tests. 


Table 5. NIST statistical test for the improved CDB-AES 


Number NIST-800-22 Tests P-value Assessment 
1. Frequency (Monobit) 0.79412 Success 
2: Block Frequency (m=128) 0.10383 Success 
3. Run 0.67902 Success 
4. Discrete Fourier Transformation 0.2932 Success 
5. Non-Overlapping Template 0.42399 Success 
6. Overlapping Template 0.6876 Success 
‘Ts Approximate Entropy 0.8941 Success 
8. Long Run 0.5113 Success 
9. Rank 0.696 Success 
10. Serial -1 0.4218 Success 

Serial -2 0.3571 i 
11. Linear Complexity 0.6955 Success 
12. Accumulative sums (forward) 0.53789 Scoss 

Accumulative sums (reverse) 0.60178 i 
13. Universal 0.83467 Success 
14. Random Excursions 0.87945 Success 
15. Random Excursions variant 0.47945 Success 


5. CONCLUSION 

This work introduced a 3D hyperchaotic system oscillation that is two-order autonomous and 
consisted of a nine-term and symmetric oscillation w.r.t x-axis used in security system with some properties 
that proved that the proposed system is hyperchaotic for getting more securing AES with a new XOR 
matrix;¢.16 (S-box) value. A novel suggestion for improving the security of the standard AES between sender 
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and receiver as a chaotic secure communication algorithm is implemented efficiently. This method is called 
CBS-AES, however, this paper generated matrix,¢,;¢ (S-box) have large keyspace and high nonlinear values. 
SAC, nonlinearity, time analysis, and NIST test were proved. The obtained matrix,¢,;¢ (S-box) demonstrated 
that all of the requirements for a decent matrix;¢.1¢ (S-box) were satisfied. We showed that the results can 
prove a way to provide high levels of security and that lead to a high level of software quality in the proposed 
system. In future work, S-box will be replaced with a new chaotic map. 
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